AN INTERACTING PARTICLE PROCESS 
RELATED TO YOUNG TABLEAUX 
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Abstract. We discuss a stochastic particle system consisting of a two-dimensional 
array of particles living in one space dimension. The stochastic evolution bears a 
certain similarity to Hammersley's process, and the particle interaction is governed 
by combinatorics of the Young tableaux. 
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1. Introduction 

The present note is related to our papers [|J, ||, [[J, and may be viewed as a 
companion to the latter paper. 

The purpose of M and |J was to construct new models of infinite-dimensional 
Markov dynamics of representation-theoretic origin. 

In || , we studied a four-parameter family of continuous time Feller Markov pro- 
cesses whose state space is the dual object U (oo) to the infinite-dimensional unitary 

group U(oo). The points of the space [/(oo) parameterize the extremal characters 
of C/(oo). Equivalently, the same space can also be viewed as the boundary dGT of 
the Gelfand-Tsetlin graph. 

In , we described a similar construction related to the infinite symmetric group 
<S(oo) that leads to a two-parameter family of Feller Markov processes on an infinite- 
dimensional cone Q C M.°° x R°° x R, called the Thoma cone. The latter space is 
an amplified version of the dual object to S'(oo): the base of the Thoma cone is the 

Thoma simplex Q = S(oo). 

l 
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The paper || explains a striking resemblance between the character theory of 
the two groups, U(oo) and S'(oo): we show how the characters of the latter group 
arise as a degeneration of those of the former group. In particular, we explain in 
what sense the Thoma cone can be viewed as a degeneration of <9GT. Moreover, the 
Markov dynamics on the Thoma cone can also be viewed as a degeneration of the 
Markov dynamics on <9GT. 

The Markov processes from and || may be viewed as interacting particle 
processes with nonlocal (or long-range) interaction. On the other hand, as shown in 
||, each of the Markov processes on <9GT is "covered" by a certain Markov process 
with local interaction, living on the path space of the graph GT. In the present 
informal note we speculate that similar covering processes should exist in the case 
of the Markov dynamics on the Thoma cone, too. 

Although we do not have a complete proof yet, we could not resist the temptation 
to introduce this model, because it provides an alternative approach to the results 
of[|. 

The note contains no proofs and presents a collection of conjectures and claims 
based on some heuristic arguments. 

2. Generalized standard Young tableaux 

Denote N := {1, 2, . . . } and N 2 := NxN. We endow N 2 with the structure of poset 
such that one element (i, j) of N 2 is declared to be greater than another element 
if i > i', j > j', and at least one of these inequalities is strict. By a (finite 
or infinite) Young diagram we mean a (finite or infinite) subset of N 2 such that if it 
contains a given element of the poset N 2 then it also contains all smaller elements. 

Let Y denote the set of all finite Young diagrams including the empty one. We 
denote a Young diagram by the letter A and write |A] for the size of A. (Infinite 
Young diagram will appear later on.) 

Recall that a standard tableau r of a given shape A G Y is a bijection r : A — > 
{1,2, ... , |A|} compatible with the partial order on A inherited from N 2 . That is, 
r(i,j) > r(i',j f ) whenever is greater than (i',f). 

Let M>o stand for the set of strictly positive real numbers. 

Definition 2.1. Following j5|, we define a generalized standard Young tableau of a 
given shape A G Y as a function h : A — > M>o such that 

• h is compatible with the partial order on A, that is, h(i, j) > h(i' , j') whenever 
(z, j) is greater than 

• the numbers h(i,j) are pairwise distinct. 

The conventional standard tableaux of a shape A G Y encode all possible ways of 
building A step by step, starting with the empty set and adding an element (i, j) at 
each step. Likewise, a generalized standard diagram encodes a similar process, only 
we imagine that elements are added at some moments hi < h% < ■ ■ • < h\\\ on the 
time axis M>o- 
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Definition 2.2. Given r G M>o and A G Y, we denote by H r (X) the set of all 
generalized standard tableaux h of the shape A, such that h(i,j) < r for all G A. 
We agree that for the empty diagram 0, the set H r (0) is a singleton. Next, we set 

H r := \jH r (\). 

AeY 

Given h G H r , we write sh(h) for the corresponding Young diagram A. It is 
convenient to define the "height" h(i,j) outside A, by setting h(i,j) = r for all 
G N 2 \ A. Then we may interpret the elements h G H r as the "height functions" 
h : N 2 — > (0,r] subject to the following conditions: 

• h(i,j) < h(i,j + 1) unless h(i,j) = h(i,j + 1) = r, 

• likewise, h(i,j) < h{i + 1, j) unless h(i,j) = h{i + 1, j) = r, 

• the set sh(h) := G N 2 : h(i,j) < r} is finite, 

• the numbers h(i,j), where ranges over sh(/i), are pairwise distinct. 

3. Markov dynamics on generalized tableaux 

As above, we fix a number r 6 E. We are going to define a continuous time 
Markov process on H r . It depends on two continuous parameters z and z' subject 
to the following condition: 

Condition 3.1. Either both parameters z and z' are complex numbers with nonzero 
imaginary part and z' = z, or both parameters are real and contained in an open 
unit interval of the form (m, m + 1) for some m G Z. 

This is equivalent to requiring that (z + k)(z' + k) > for every k G Z. In 
particular, Condition [3J] implies that 

q(i,j) := ( z + j-i)( z ' + j-i) > , V(z,j) G N 2 , (3.1) 

which is used in the sequel. 

We need one more notation. Given h G H r and (z, j) G N 2 , we set 

max{h(i — 1, j), /i(^,J — 1)} , if i > 1 and j > 1; 

/i(2-l,l), ifz>l,j = l; 

,0, ifz=j = l. 

Let us fix a couple (z, z') of parameters satisfying Condition |3.1| . We consider an 
infinite collection of independent Poisson processes indexed by the elements G 
N 2 , where each process is defined in the horizontal strip 

D r := {(t, x) G 1R 2 , -oo < t < +oo, < x < r}, (3.3) 
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and the (z,j)th process has constant rate equal to q(i,j) (see (|3.1|) ). We denote 
by n r (i,j) C D r the random point configuration from the (i,j)th process. We may 
assume that n r (i,j) has at most one point on each vertical line. 

Definition 3.2. We introduce a Markov process Xr ' on H r , whose time evolution 
is composed of two components, one deterministic and the other one stochastic: 

• The deterministic component prescribes the height h(i,j) (for every (i,j) £ 
sh(/i)) to grow with the varying velocity v(y) := y(y + 1) depending on the 
current value y := h(i,j), until the height attains the maximum value r, 
where it stabilizes. 

• The stochastic component, on the contrary, forces the height to instantly 
drop to some level. More precisely, given £ N 2 , the height h(i,j) drops 
precisely at those time moments t for which the vertical line with abscissa t 
contains a point (t,x) £ n r (i,j) C D r such that the ordinate x satisfies the 
constraints h}(i,j) < x < h(i,j). Then the height instantly takes the new 
value equal to x. 

Note that the deterministic evolution preserves the set H r : this follows from the 
very definition of H r and the fact that the velocity v(x) is a strictly increasing 
function in x. Likewise, the stochastic evolution also preserves H r , by virtue of the 
requirement that for the drop to take place x must be greater than h^(i,j). More 
pedantically, we had to include in the above definition the additional requirement 
that a jump h(i,j) — > x is forbidden whenever x is equal to h(i',f) for some £ 
sh(/i) distinct from but such an event has probability 0. 

We may interpret Xr ' as an evolution of a system of particles on (0, r] indexed 
by the two-dimensional array N 2 . It is worth noting that the stochastic component of 
the process includes an interaction between the particles, because the lower bound 
h^(i,j) for the possible drop of the (i,j)th particle depends on the positions of 
the neighboring particles. Here the "neighbors" are defined as the particles with 
neighboring indices (i — l,j) and (i,j — 1). This kind of interaction is similar to 
that encountered in the exclusion processes. 

On the other hand, because of the presence of two components, deterministic and 
stochastic, the process Xr may be viewed clS db member of the class of piecewise 
deterministic Markov processes introduced by Davis J7|. 

The following claim looks very plausible to us. 

Conjecture 3.3. The process x\ z,z ^ on H r does not explode and so has infinite life 
time almost surely. 

4. Evolution of Gibbs measures 

Here we discuss a connection between the process Xr Z ' z ^ just defined and the 
process Xr Z,z ^ discussed in our papers and 0. 
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Let us recall the definition of Xr Z,z ^ (see || Definition 8.9]). This is a continuous 
time jump Markov process with the state space Y. The evolution is given the 
following Q-matrix (the matrix of jump rates): 

Qi'-'H\, A + □) = r(z + c(D))(z' + C (D))^^±EL, □ e \\ 

^>(A,A-D) = (r+ l) |A|d 7 ( V D) . O.A-, ^ 

dim A 

-Qi z ' z '\X,X) = (2r + l)\\\+rzz'. 

The notation is the following: A is a Young diagram; A ± □ is another diagram 
obtained from A by appending/removing a box □; A 1 * 1 is the set of those boxes that 
can be appended to (respectively, removed from) A; Q(A, A ± □) is the rate of the 
jump A — > A ± □; finally, — Q(X, A) is equal to the sum of the rates of all possible 
jumps A — > A ± □. 

As shown in ||, this Q- matrix is regular, meaning that the corresponding jump 
Markov process does not explode. 

Let us return to the sets H r (X) introduced in Definition |2.2| . For every A G Y, let 
Tab(A) denote the finite set consisting of all (conventional) standard tableaux of the 
shape A. As seen from Definition |3.2| , there is a natural projection H r (X) — > Tab(A). 
For every standard tableau r G Tab(A), the fiber of this projection over r can be 
viewed as an open simplex A(r) of dimension N = |A|, formed by the ordered N- 
tuples of reals (xi, . . . , xjv) such that < x\ < ■ ■ ■ < xn < r. It follows that H r (X) 
can be viewed as a bounded open subset of M. N whose closure H r (X) is a convex 
polytope endowed with a triangulation. 

For instance, if A = (2, 1), then N = 3, H r (X) consists of the triples (xi, x 2 , x 3 ) = 
(/i(l, 1), h(l, 2), h(2, 1)) subject to the conditions 

< Xi < r, < x 2 < r, < x 3 < r, Xi < x 2 , X\< x 3 , x 2 ^ £ 3 , 



and H r (X) is obtained by removing the last inequality and making the remaining 
inequalities weak. The polytope H r (X) is the union of two closed simplices, which 
are singled out by the inequalities x 2 < x% and x 2 > x^, respectively. 

Definition 4.1. Recall (see Definition |2.2|) that H r is the disjoint union of the sets 
H r (X), where A ranges over Y. Following we define a Gibbs measure on H r as a 
probability measure such that its restriction to each subset H r (X) is proportional to 
the Lebesgue measure. The set of all Gibbs measures is denoted by 

Obviously, the natural projection H r — > Y establishes a one-to-one correspondence 
between 5f r and the set ^#(Y) of all probability measures on Y. 

Claim 4.2. The process Xr ' preserves the set & r of Gibbs measures, and the 
evolution of the Gibbs measures induced by the process Xr coincides, under the 
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bijection & r — > ^%(Y), with the evolution of the probability measures on Y induced 
by the process X^ z ' z ^ . 

Actually, we can rigorously prove only the infinitesimal version of the claim: the 
application of the infinitesimal generator of Xr ' to a Gibbs measure translates to 
the application of the generator of Xr Z ' z to the corresponding measure on Y. 

Remark 4.3. Here is a simple yet curious formal identity used in the proof. Assume 
we are given a standard tableau of the shape N 2 , that is, a total order on the set N 2 
compatible with its partial order. Next, for every G N 2 , set 



'max{(i - l,j), - 1)} , if % > 1 and j > 1, 

(U-i), if i = l, j > l, 

(2-1,1), ifi>l, j = l, 

undefined, if = (1,1), 



and 



(h j) f := min{(i + l,j), + 1)} , 

where the maximum and minimum are taken relative to the prescribed total order 
on N 2 . Finally, attach to every G N 2 a formal variable y(i,j). Then the identity 
in question is 

(y((i,j) t ) + y((i,3) i )-2y(i,j))(z + j-i)(z' + j-i) = 2 y&j) 

(i,j)GN 2 («J)eN 2 
with the agreement that y((l, 1)^) := 0. 

rn (z z') 

As explained in |6[ , the process X. r ' has a stationary distribution, the so called 
mixed z-measure, denoted by Mr Z ' z \ Let Mr Z ' z ^ denote the corresponding Gibbs 
measure on H r . Claim |L2| implies that M^ z ' z>) serves as the stationary distribution 
for the process Xr Z ' z . 

Claim says that the process Xr Z ' z ' in some sense "covers" the jump process 

(z z f ) 

Xr . Note that Young diagrams A G Y can be represented as particle configura- 
tions, so that Xr Z ' z \ like X r z ' z \ can also be interpreted as an interacting particle 
process. However, a substantial difference between Xr Z ' z ^ and Xj- Z ' z ' is that the 
particle interaction is local in the former process and highly non-local in the latter 
one. In this sense, Xr Z ' z ' seems to be simpler than X^ z,z \ 

We showed in [BJ that for every pair r' > r there exists a "link" K r r (an infinite 
stochastic matrix of format Yx Y), which intertwines the processes X { r f' z,) and X^ z ' ] 
for all (z, z'). For the covering processes the picture is simpler in the sense that the 
processes X {z r z,) and X^ z ' ] are "linked" by an ordinary map, the truncation map 
( |5.1| ) that we now define. 
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5. Dynamics on infinite generalized tableaux 

Given r' > r, we define the truncation map H r r — > H r as the transform 

h{i,j) -+mm(h(i,j),r), V(iJ) £ N 2 . (5.1) 

Obviously, for a triple r" > r' > r, the composed map H r n — > H r i — > H r is the same 
as H r ii — > H r . Therefore, we may define the projection limit space 

H:=\^mH r , r — > +oo. 

The elements of H can be viewed as the functions h : N 2 — > M>o U {+00} such that 

• h(i,j) < h{i + 1, j) unless h(i,j) = h{i + 1, j) = +00; 

• h(i,j) < h{i,j + 1) unless h(i,j) = h(i,j + 1) = +00; 

• the finite values h(i,j) are pairwise distinct and do not have accumulation 
points on M>o- 

Note that the natural projection H — > H r is still given by (|5.1|). Note also that, 
for h G H, the set 

sh(h) := e N 2 : h(i, j) < +00} C N 2 

is a (possibly infinite) Young diagram. 

Every element h £ H may be interpreted as a system of particles on the extended 
halfline R >0 U {+00}, indexed by N 2 (we assume that the (i,j)th particle has coor- 
dinate x = h(i,j)). If sh(h) = N 2 , then all particles are on K>o, but it may happen 
that sh(h) is a proper subset of N 2 ; then there is an infinite reservoir of particles at 
infinity, indexed by the elements of N 2 \ sh(h). An important requirement included 
in the definition is that the particle configuration on IR>o is always locally finite 
meaning that there are finitely many particles in any bounded interval. 

Claim 5.1. The family of Markov processes {X { r Z ' Z,) : r £ M>o} is consistent with 
the truncation maps H r > — > H r and so there exists a unique Markov process X^- Z ' z '^ on 
the space H , consistent with this family by means of the truncation maps H — >■ H r . 

In other words, X^ z,z ) can be viewed as the projective limit of the processes Xr ' ^ . 

The Markov dynamics of such a system, given by the process X^ z,z '\ is described 
in the same way as for the truncated processes Xr z ' z , see Section ^| above, only the 
open interval (0,r) should be replaced by the halfline M>o, and the endpoint r is 
shifted to +00. 

Informally, in the particle system interpretation, the deterministic component 
looks as the accelerating movement of particles to the right with velocity v(y) = y(y+ 
1), while the the concurrent stochastic component forces the particles to instantly 
jump to the left. The latter component is driven by a two-dimensional array {7r(z, j) : 
£ N 2 } of independent Poisson processes in the half-plane D := {(t,x) : t £ 
M, x > 0}, where the (i,j)th process has constant rate q(i,j) given by (|3.1|) . Note 
that the particles can escape from IR >0 to infinity and return back to M>o- 
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Note also that the sample trajectories of X^ z ' z "> look more sophisticated than 
those of the truncated process X^ z ' z \ This is due to a major difference between 
the Poisson processes in the strip D r and in the half-plane D: the pushforward of a 
Poisson configuration under the projection to the i-axis is locally finite in the former 
case but not in the latter case. 

Now we extend our definition of Gibbs measures to the space H. 

Definition 5.2. A probability measure on H is said to be a Gibbs measure if its 
pushforward under the truncation map H — > H r is a Gibbs measure on H r in the 



sense of Definition |4.1| , for every r G K>o- The set of all Gibbs measures on H is 
denoted by 3f . 

The above discussion shows that the process X°° preserves the Gibbs measures. 

On the other hand, as explained in ||, there is a one-to-one correspondence 
£f <H- ^(Q), where Q is an infinite-dimensional cone in R°° x K°° x R, called the 
Thoma cone, and ^#(0) denotes the space of probability measures on Q. Further, 
we showed in || that the family {Xr Z ' z : r G R>o} of Markov processes on Y 
determines, via the links A£ , a Markov process X^ z ' z "> on the Thoma cone Vt. Because 
of the bijection -H- ^#(0), this leads to the following conclusion: 

Claim 5.3. The process X^ z,z > covers the process X^ in the sense that the evo- 
lution of arbitrary probability measures on the Thoma cone induced by X^ z,z > is the 
same as the evolution of the Gibbs measures on H induced by X^ z ' z '\ 

Here is a reformulation. Given a point u G O, let us denote by the Gibbs 
measure on H corresponding to the Dirac measure at u. The correspondence u i— > 
G^ can be viewed as a Markov kernel = A^(w, dh). Then the above claim means 
that A| intertwines X^ z ^ and X^ z ' z '\ 

In 0, we showed that the process A* 2 ' 2 ') on the Thoma cone has a stationary 
distribution M^ z,z '\ which is defined by the family {Mr Z ' z : r G M>o} via the links 
A£ . Let G^ z ' z ' stand for the corresponding Gibbs measure on H; in other words, 
q{z,z') j g transform of M^ z,z '^ by A^. 

Claim 5.4. The Gibbs measure serves as a stationary distribution for X^ z,z '. 

6. Concluding remarks 

Besides the truncation maps H — > H r there exist other projections under which 
the Markov property of the process X^ z ' z '^ is not destroyed. Namely, one may fix 
an arbitrary Young diagram, finite or infinite, and focus on those particles that are 
indexed by the elements of that diagram. 

The simplest such example is obtained when the diagram in question consists of 
the single element (1,1). Then the model represents a single particle that moves on 
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the extended halfline M>o U {+00}. The dynamics is described exactly as before: 
a deterministic movement to the right with the velocity v(y) = y(y + 1) combined 
with instant jumps to the left directed by the Poisson process in the half-plane D 
with constant rate q(l, 1) = zz' . 

A more complicated model arises when we take the infinite one-row diagram 
'■ j G N}. Then we have countably many particles, enumerated by the 
numbers j 6 N and occupying positions x\ < X2 < ■ ■ ■ on M >0 with a possible 
infinite reservoir at infinity, meaning that it may happen that Xj = +00 for all 
j large enough. The dynamics in this model has a certain resemblance with the 
Hammersley process studied in Aldous-Diaconis |JT| . 

Finally, note that the process X^ z ' z '^ admits a Plancherel-type degeneration in the 
spirit of H Section 10]. The idea is to scale the space variable x by substituting 
x{zz')~ l instead of x, and pass to a limit as both z and z' go to +00. Then we get a 
simplified model, where the velocity in the deterministic component depends on the 
coordinate y G M>o via v(y) = y, and the backward jump rules are the same for all 
particles (they are governed by independent Poisson processes in D with constant 
rate q = 1). 
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